Refining the Global Phylogeny of Mitochondrial N1a, X, and HV2 Haplogroups Based on Rare Mitogenomes from Croatian Isolates.

Mitochondrial DNA (mtDNA) has been used for decades as a predominant tool in population genetics and as a valuable addition to forensic genetic research, owing to its unique maternal inheritance pattern that enables the tracing of individuals along the maternal lineage across numerous generations. The dynamic interplay between evolutionary forces, primarily genetic drift, bottlenecks, and the founder effect, can exert significant influence on genetic profiles. Consequently, the Adriatic islands have accumulated a subset of lineages that exhibits remarkable absence or rarity within other European populations. This distinctive genetic composition underscores the islands' potential as a significant resource in phylogenetic research, with implications reaching beyond regional boundaries to contribute to a global understanding. In the initial attempt to expand the mitochondrial forensic database of the Croatian population with haplotypes from small isolated communities, we sequenced mitogenomes of rare haplogroups from different Croatian island and mainland populations using next-generation sequencing (NGS). In the next step and based on the obtained results, we refined the global phylogeny of haplogroup N1a, HV2, and X by analyzing rare haplotypes, which are absent from the current phylogenetic tree. The trees were based on 16 novel and 52 previously published samples, revealing completely novel branches in the X and HV2 haplogroups and a new European cluster in the ancestral N1a variant, previously believed to be an exclusively African-Asian haplogroup. The research emphasizes the importance of investigating geographically isolated populations and their unique characteristics within a global context.

Stable methylation loci are associated with systolic blood pressure in a Croatian island population.

Background: The objective was to identify stable and dynamic DNA methylation loci associated with cardiometabolic traits among an adult population from the Croatian island of Hvar. Materials & methods: An epigenome-wide association study was conducted using peripheral blood longitudinally collected at two time points 10 years apart via Infinium MethylationEPIC beadarray (n = 112). Stable and dynamic loci were identified using linear mixed models. Associations between cardiometabolic traits and loci were assessed using linear models. Results: 22 CpG loci were significantly associated with systolic blood pressure. Twenty were stable and two were dynamic. Conclusion: Multiple genes may be involved in the determination of systolic blood pressure level via stable epigenetic programming, potentially established earlier in life.

Cardiovascular disease is the leading cause of death worldwide. Previous studies have found that genetics incompletely explain susceptibility to cardiovascular disease. To find new potential risk factors, the authors investigated the possible contribution of DNA methylation (modifications to DNA that can affect gene expression but do not alter the underlying genetic code) in an adult population on the Croatian island of Hvar, which has a high number of people with cardiovascular and metabolic disease. By examining DNA methylation in blood collected at two time points, 10 years apart, the authors were able to identify DNA methylation that either stayed the same over time (stable) or changed the most over time (dynamic). These were then compared with clinical test results related to cardiovascular or metabolic diseases to determine if they are associated. Twenty-two methylation sites were found to be associated with systolic blood pressure. Of those, 20 were considered stable and two were dynamic. Additionally, there was one stable methylation site associated with serum calcium and one with C-reactive protein. These findings suggest that systolic blood pressure may be regulated through stable DNA methylation that is potentially established earlier in life.

Genetic sub-structuring of Croatian island populations in the Southeastern European context: a meta-analysis.

AIM: To use the method of meta-analysis to assess the influence of island population isolation on the sub-structuring of the Croatian population, as well as the influence of regional population groups on the sub-structuring of the Southeastern European population with regard to basic population genetic statistical parameters calculated by using STR locus analysis. METHODS: Bio-statistical analyses were performed for 2877 unrelated participants of both sexes from Southeastern Europe. Nine autosomal STR loci (D3S1358, vWA, FGA, TH01, TPOX, CSF1PO, D5S818, D13S317, and D7S82) were analyzed by using standard F-statistics and population structure analysis (Structure software). RESULTS: Genetic differentiation of Croatian subpopulations assessed with the FST method was higher at the level of the Croatian population (0.005) than at the level of Southeastern Europe (0.002). The island of Vis showed the most pronounced separation in the Croatian population, and Albanians from Kosovo in the population of Southeast Europe, followed by Croatia, Bosnia and Herzegovina, and Hungary. CONCLUSION: The higher structure of Croatian subpopulations in relation to Southeastern Europe suggest a certain degree of genetic isolation, most likely due to the influence of endogamy within rural island populations.

Influence of genetic substructuring of statistical forensic parameters on genetic short tandem repeat markers in the populations of Southeastern Europe.

AIM: To investigate the influence of specific intrapopulation genetic structures on interpopulation relationships. Special focus was the influence of island population isolation on the substructuring of the Croatian population, and the influence of regional population groups on the substructuring of Southeast European populations. METHODS: Autosomal short tandem repeat (STR) loci were analyzed by using four forensic parameters: matching probability (PM), power of discrimination (PD), power of exclusion (PE), and polymorphic information content (PIC) on a sample of 2877 unrelated participants of both sexes. A sample set comprising 590 participants was analyzed for the first time, and 2287 participants were included from previous studies. The analysis was performed with PowerStats v. 1.2. RESULTS: The analysis of forensic parameters for all nine loci in the Croatian subpopulations showed the largest deviations in the populations of the islands of Korcula and Hvar. The smallest deviations were found in the mainland population. As for Southeast European populations, the largest deviations were found in the population of North Macedonia, followed by Romania, Albanians from Kosovo, and Montenegro, while the smallest deviations were found in the population of Hungary. CONCLUSION: The comparison of forensic parameters between different subpopulations of Croatia and Southeast Europe indicates that the isolation of individual Croatian subpopulations and rare alleles in their gene pool affect the values of forensic parameters. Specific features of (sub)populations should be taken into account for appropriate sampling of the total population when creating a DNA database of STR markers.

Genetic landscape of the Hvar Island - highlight of a 50-year long bioanthropological research.

The paper presents an overview of the 50-year long bioanthropological research of the Hvar islanders and depicts the maternal and paternal genetic landscape of the Hvar population (mtDNA and NRY lineages) in more detail. MtDNA haplogroups were determined in 169 and NRY haplogroups in 407 autochthonous individuals from the Hvar Island. The relatively high level of diversity of mtDNA and NRY lineages has been observed, however with interesting deviations from both the maternal (F1b1 lineage) and paternal (Q2a1a lineage) perspective. Additionally, population substructuring revealed differences between Hvar communities (east-west substructuring), in line with the ethnohistoric background and observed migration patterns on the island. Genetic analysis of the Hvar islanders presents a highlight of the 50-year long anthropological research on this island and offers insight into the current genetic structure of Dalmatia, Croatia, shaped by dynamic and diverse population movements throughout history.

Association of perfluoroalkyl substances exposure with cardiometabolic traits in an island population of the eastern Adriatic coast of Croatia.

BACKGROUND: Exposure to perfluoroalkyl substances (PFAS), ubiquitous environmental contaminants, may be related to cardiometabolic diseases in adults. Studies in European populations to examine the association of PFAS exposure and comprehensive cardiometabolic traits and metabolic syndrome (MetS) are limited. METHODS: In this pilot cross-sectional study of a well-characterized adult population of the island of Hvar, situated off the eastern Adriatic coast of Croatia, we measured PFAS concentrations in plasma samples collected during 2007-2008 and examined their cross-sectional associations with cardiometabolic traits and MetS after adjustment of covariates (n = 122). PFAS investigated in this study included perfluorooctane sulfonic acid (PFOS), perfluorooctanoic acid (PFOA), perfluorohexane sulfonic acid (PFHxS), and perfluorononanoic acid (PFNA). RESULTS: The geometric mean (range) was 8.91 (2.36, 33.67) ng/mL for PFOS, 2.87 (1.03, 8.02) ng/mL for PFOA, 0.77 (0.25, 2.40) ng/mL for PFHxS, and 1.29 (0.48, 3.46) ng/mL for PFNA, with frequency of detection at 100%, 100%, 95.9%, and 100%, respectively. PFOS, PFOA, and PFNA concentrations were positively associated with the risk of MetS as defined by the Adult Treatment Panel III (ATP III) criteria, with estimated odds ratios and 95% confidence intervals at 1.89 (0.93, 3.86), 2.19 (0.88, 5.44), and 2.95 (1.12, 7.80), respectively, with only PFNA reaching statistical significance. PFNA concentrations were associated with increased risk of overweight or obesity. CONCLUSIONS: Background exposure to PFOS, PFOA, and PFNA was marginally associated with increased risk of MetS in this small study, and these results should be confirmed with a larger sample size and longitudinal follow-up.

Reconstructing the genetic history of late Neanderthals.

Although it has previously been shown that Neanderthals contributed DNA to modern humans, not much is known about the genetic diversity of Neanderthals or the relationship between late Neanderthal populations at the time at which their last interactions with early modern humans occurred and before they eventually disappeared. Our ability to retrieve DNA from a larger number of Neanderthal individuals has been limited by poor preservation of endogenous DNA and contamination of Neanderthal skeletal remains by large amounts of microbial and present-day human DNA. Here we use hypochlorite treatment of as little as 9 mg of bone or tooth powder to generate between 1- and 2.7-fold genomic coverage of five Neanderthals who lived around 39,000 to 47,000 years ago (that is, late Neanderthals), thereby doubling the number of Neanderthals for which genome sequences are available. Genetic similarity among late Neanderthals is well predicted by their geographical location, and comparison to the genome of an older Neanderthal from the Caucasus indicates that a population turnover is likely to have occurred, either in the Caucasus or throughout Europe, towards the end of Neanderthal history. We find that the bulk of Neanderthal gene flow into early modern humans originated from one or more source populations that diverged from the Neanderthals that were studied here at least 70,000 years ago, but after they split from a previously sequenced Neanderthal from Siberia around 150,000 years ago. Although four of the Neanderthals studied here post-date the putative arrival of early modern humans into Europe, we do not detect any recent gene flow from early modern humans in their ancestry.

A high-coverage Neandertal genome from Vindija Cave in Croatia.

To date, the only Neandertal genome that has been sequenced to high quality is from an individual found in Southern Siberia. We sequenced the genome of a female Neandertal from ~50,000 years ago from Vindija Cave, Croatia, to ~30-fold genomic coverage. She carried 1.6 differences per 10,000 base pairs between the two copies of her genome, fewer than present-day humans, suggesting that Neandertal populations were of small size. Our analyses indicate that she was more closely related to the Neandertals that mixed with the ancestors of present-day humans living outside of sub-Saharan Africa than the previously sequenced Neandertal from Siberia, allowing 10 to 20% more Neandertal DNA to be identified in present-day humans, including variants involved in low-density lipoprotein cholesterol concentrations, schizophrenia, and other diseases.

Neandertal and Denisovan DNA from Pleistocene sediments.

Although a rich record of Pleistocene human-associated archaeological assemblages exists, the scarcity of hominin fossils often impedes the understanding of which hominins occupied a site. Using targeted enrichment of mitochondrial DNA, we show that cave sediments represent a rich source of ancient mammalian DNA that often includes traces of hominin DNA, even at sites and in layers where no hominin remains have been discovered. By automation-assisted screening of numerous sediment samples, we detected Neandertal DNA in eight archaeological layers from four caves in Eurasia. In Denisova Cave, we retrieved Denisovan DNA in a Middle Pleistocene layer near the bottom of the stratigraphy. Our work opens the possibility of detecting the presence of hominin groups at sites and in areas where no skeletal remains are found.

Genetic heritage of Croatians in the Southeastern European gene pool-Y chromosome analysis of the Croatian continental and Island population.

OBJECTIVES: The research objective of this study is to enlarge and deepen the Y chromosome research on the Croatian population and enable additional insights into the population diversity and historic events that shaped the current genetic landscape of Croatia and Southeastern Europe (SEE). MATERIALS AND METHODS: A high-resolution phylogenetic and phylogeographic analysis of 66 biallelic (SNPs) and 17 microsatellite (STRs) markers of the Y chromosome was performed using 720 Croatian samples. The obtained results were placed in a wider European context by comparison with ∼4450 samples from a number of other European populations. RESULTS: A high diversity of haplogroups was observed in the overall Croatian sample, and all typical European Y chromosome haplogroups with corresponding clinal patterns were observed. Three distinct genetic signals were identifiable in the Croatian paternal gene pool - I2a1b-M423, R1a1a1b1a*-M558, and E1b1b1a1b1a-V13 haplogroups. DISCUSSION: The analyses of the dominant and autochthonous I2a1b-M423 lineage (>30%) suggest that SEE had a significant role in the Upper Paleolithic, the R1a1a1b1a*-M558 lineage (19%) represents a signal from present day Slavic populations of Central Europe in the Croatian population, and the phylogeography of the E1b1b1a1b1a-V13 clade (around 9%) implies cultural diffusion of agriculture into Europe via the Balkan Peninsula. Am. J. Hum. Biol., 2016. © 2016 Wiley Periodicals, Inc. Am. J. Hum. Biol. 28:837-845, 2016. © 2016Wiley Periodicals, Inc.

Ancient gene flow from early modern humans into Eastern Neanderthals.

It has been shown that Neanderthals contributed genetically to modern humans outside Africa 47,000-65,000 years ago. Here we analyse the genomes of a Neanderthal and a Denisovan from the Altai Mountains in Siberia together with the sequences of chromosome 21 of two Neanderthals from Spain and Croatia. We find that a population that diverged early from other modern humans in Africa contributed genetically to the ancestors of Neanderthals from the Altai Mountains roughly 100,000 years ago. By contrast, we do not detect such a genetic contribution in the Denisovan or the two European Neanderthals. We conclude that in addition to later interbreeding events, the ancestors of Neanderthals from the Altai Mountains and early modern humans met and interbred, possibly in the Near East, many thousands of years earlier than previously thought.

Variability and distribution of COL1A2 (type I collagen) polymorphisms in the central-eastern Mediterranean Basin.

BACKGROUND: The most abundant of the collagen protein family, type I collagen is encoded by the COL1A2 gene. The COL1A2 restriction fragment length polymorphisms (RFLPs) EcoRI, RsaI and MspI in samples from several different central-eastern Mediterranean populations were analysed and found to be potentially informative anthropogenetic markers. AIM: The objective was to define the genetic variability of COL1A2 in the central-eastern Mediterranean and to shed light on its genetic distribution in human groups over a wide geographic area. SUBJECTS AND METHODS: PCR-RFLP analysis of EcoRI, RsaI and MspI polymorphisms of the COL1A2 gene was performed on oral swab and blood samples from 308 individuals from the central-eastern Mediterranean Basin. The genetic similarities among these groups and other populations described in the literature were investigated through correspondence analysis. RESULTS: Single-marker data and haplotype frequencies seemed to suggest a genetic homogeneity within the European populations, whereas a certain degree of differentiation was noted for the Egyptians and the Turks. CONCLUSIONS: The genetic variability in the central-eastern Mediterranean area is probably a result of the geographical barrier of the Mediterranean Sea, which separated European and African populations over time.

Genetic Heritage of the Balto-Slavic Speaking Populations: A Synthesis of Autosomal, Mitochondrial and Y-Chromosomal Data.

The Slavic branch of the Balto-Slavic sub-family of Indo-European languages underwent rapid divergence as a result of the spatial expansion of its speakers from Central-East Europe, in early medieval times. This expansion-mainly to East Europe and the northern Balkans-resulted in the incorporation of genetic components from numerous autochthonous populations into the Slavic gene pools. Here, we characterize genetic variation in all extant ethnic groups speaking Balto-Slavic languages by analyzing mitochondrial DNA (n = 6,876), Y-chromosomes (n = 6,079) and genome-wide SNP profiles (n = 296), within the context of other European populations. We also reassess the phylogeny of Slavic languages within the Balto-Slavic branch of Indo-European. We find that genetic distances among Balto-Slavic populations, based on autosomal and Y-chromosomal loci, show a high correlation (0.9) both with each other and with geography, but a slightly lower correlation (0.7) with mitochondrial DNA and linguistic affiliation. The data suggest that genetic diversity of the present-day Slavs was predominantly shaped in situ, and we detect two different substrata: 'central-east European' for West and East Slavs, and 'south-east European' for South Slavs. A pattern of distribution of segments identical by descent between groups of East-West and South Slavs suggests shared ancestry or a modest gene flow between those two groups, which might derive from the historic spread of Slavic people.

Genetic structure and admixture between Bayash Roma from northwestern Croatia and general Croatian population: evidence from Bayesian clustering analysis.

The European Roma represent a transnational mosaic of minority population groups with different migration histories and contrasting experiences in their interactions with majority populations across the European continent. Although historical genetic contributions of European lineages to the Roma pool were investigated before, the extent of contemporary genetic admixture between Bayash Roma and non-Romani majority population remains elusive. The aim of this study was to assess the genetic structure of the Bayash Roma population from northwestern Croatia and the general Croatian population and to investigate the extent of admixture between them. A set of genetic data from two original studies (100 Bayash Roma from northwestern Croatia and 195 individuals from the general Croatian population) was analyzed by Bayesian clustering implemented in STRUCTURE software. By re-analyzing published data we intended to focus for the first time on genetic differentiation and structure and in doing so we clearly pointed to the importance of considering social phenomena in understanding genetic structuring. Our results demonstrated that two population clusters best explain the genetic structure, which is consistent with social exclusion of Roma and the demographic history of Bayash Roma who have settled in NW Croatia only about 150 years ago and mostly applied rules of endogamy. The presence of admixture was revealed, while the percentage of non-Croatian individuals in general Croatian population was approximately twofold higher than the percentage of non-Romani individuals in Roma population corroborating the presence of ethnomimicry in Roma.

The phylogenetic and geographic structure of Y-chromosome haplogroup R1a.

R1a-M420 is one of the most widely spread Y-chromosome haplogroups; however, its substructure within Europe and Asia has remained poorly characterized. Using a panel of 16 244 male subjects from 126 populations sampled across Eurasia, we identified 2923 R1a-M420 Y-chromosomes and analyzed them to a highly granular phylogeographic resolution. Whole Y-chromosome sequence analysis of eight R1a and five R1b individuals suggests a divergence time of ∼25,000 (95% CI: 21,300-29,000) years ago and a coalescence time within R1a-M417 of ∼5800 (95% CI: 4800-6800) years. The spatial frequency distributions of R1a sub-haplogroups conclusively indicate two major groups, one found primarily in Europe and the other confined to Central and South Asia. Beyond the major European versus Asian dichotomy, we describe several younger sub-haplogroups. Based on spatial distributions and diversity patterns within the R1a-M420 clade, particularly rare basal branches detected primarily within Iran and eastern Turkey, we conclude that the initial episodes of haplogroup R1a diversification likely occurred in the vicinity of present-day Iran.

Fish and Shellfish Intake and Diabetes in a Costal Population of the Adriatic.

The objective of the study was to examine the association between fish and shellfish intake and diabetes in an island population, and the design of the study was Cross-sectional. Two independent population-based field surveys were conducted in Hvar Island of the eastern Adriatic coast of Croatia in May 2007 and May 2008, with a total of 1,379 adult participants. In multivariable logistic regression models, total fish intake was positively associated with diabetes prevalence in the total population (OR(Q4 vs. Q1) = 1.64; 95% CI = 1.01-2.66; p-trend = 0.09). Oily fish intake also exhibited a positive association with diabetes prevalence in the total population (OR(Q4 vs. Q1) = 2.22; 95% CI = 1.35-3.64; p-trend = 0.01) and in analyses stratified by body mass index, males and those with a high waist circumference. The study suggests an association between oily fish intake and diabetes in the population of the Hvar Island in Croatia. Longitudinal studies incorporating measures of persistent organic pollutants and local cooking practices are warranted to identify factors in fatty fish that may influence the development or persistence of diabetes.

Standing at the gateway to Europe--the genetic structure of Western balkan populations based on autosomal and haploid markers.

Contemporary inhabitants of the Balkan Peninsula belong to several ethnic groups of diverse cultural background. In this study, three ethnic groups from Bosnia and Herzegovina - Bosniacs, Bosnian Croats and Bosnian Serbs - as well as the populations of Serbians, Croatians, Macedonians from the former Yugoslav Republic of Macedonia, Montenegrins and Kosovars have been characterized for the genetic variation of 660 000 genome-wide autosomal single nucleotide polymorphisms and for haploid markers. New autosomal data of the 70 individuals together with previously published data of 20 individuals from the populations of the Western Balkan region in a context of 695 samples of global range have been analysed. Comparison of the variation data of autosomal and haploid lineages of the studied Western Balkan populations reveals a concordance of the data in both sets and the genetic uniformity of the studied populations, especially of Western South-Slavic speakers. The genetic variation of Western Balkan populations reveals the continuity between the Middle East and Europe via the Balkan region and supports the scenario that one of the major routes of ancient gene flows and admixture went through the Balkan Peninsula.

Patterns of coding variation in the complete exomes of three Neandertals.

We present the DNA sequence of 17,367 protein-coding genes in two Neandertals from Spain and Croatia and analyze them together with the genome sequence recently determined from a Neandertal from southern Siberia. Comparisons with present-day humans from Africa, Europe, and Asia reveal that genetic diversity among Neandertals was remarkably low, and that they carried a higher proportion of amino acid-changing (nonsynonymous) alleles inferred to alter protein structure or function than present-day humans. Thus, Neandertals across Eurasia had a smaller long-term effective population than present-day humans. We also identify amino acid substitutions in Neandertals and present-day humans that may underlie phenotypic differences between the two groups. We find that genes involved in skeletal morphology have changed more in the lineage leading to Neandertals than in the ancestral lineage common to archaic and modern humans, whereas genes involved in behavior and pigmentation have changed more on the modern human lineage.

Maternal genetic heritage of Southeastern Europe reveals a new Croatian isolate and a novel, local sub-branching in the x2 haplogroup.

High mtDNA variation in Southeastern Europe (SEE) is a reflection of the turbulent and complex demographic history of this area, influenced by gene flow from various parts of Eurasia and a long history of intermixing. Our results of 1035 samples (488 from Croatia, 239 from Bosnia and 130 from Herzegovina, reported earlier, and 97 Slovenians and 81 individuals from Zumberak, reported here for the first time) show that the SEE maternal genetic diversity fits within a broader European maternal genetic landscape. The study also shows that the population of Zumberak, located in the continental part of Croatia, developed some unique mtDNA haplotypes and elevated haplogroup frequencies due to distinctive demographic history and can be considered a moderate genetic isolate. We also report seven samples from the Bosnian population and one Herzegovinian sample designated as X2* individuals that could not be assigned to any of its sublineages (X2a'o) according to the existing X2 phylogeny. In an attempt to clarify the phylogeny of our X2 samples, their mitochondrial DNA has been completely sequenced. We suppose that these lineages are signs of local microdifferentiation processes that occurred in the recent demographic past in this area and could possibly be marked as SEE-specific X2 sublineages.

Glycans are a novel biomarker of chronological and biological ages.

Fine structural details of glycans attached to the conserved N-glycosylation site significantly not only affect function of individual immunoglobulin G (IgG) molecules but also mediate inflammation at the systemic level. By analyzing IgG glycosylation in 5,117 individuals from four European populations, we have revealed very complex patterns of changes in IgG glycosylation with age. Several IgG glycans (including FA2B, FA2G2, and FA2BG2) changed considerably with age and the combination of these three glycans can explain up to 58% of variance in chronological age, significantly more than other markers of biological age like telomere lengths. The remaining variance in these glycans strongly correlated with physiological parameters associated with biological age. Thus, IgG glycosylation appears to be closely linked with both chronological and biological ages. Considering the important role of IgG glycans in inflammation, and because the observed changes with age promote inflammation, changes in IgG glycosylation also seem to represent a factor contributing to aging. SIGNIFICANCE STATEMENT: Glycosylation is the key posttranslational mechanism that regulates function of immunoglobulins, with multiple systemic repercussions to the immune system. Our study of IgG glycosylation in 5,117 individuals from four European populations has revealed very extensive and complex changes in IgG glycosylation with age. The combined index composed of only three glycans explained up to 58% of variance in age, considerably more than other biomarkers of age like telomere lengths. The remaining variance in these glycans strongly correlated with physiological parameters associated with biological age; thus, IgG glycosylation appears to be closely linked with both chronological and biological ages. The ability to measure human biological aging using molecular profiling has practical applications for diverse fields such as disease prevention and treatment, or forensics.